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Asteroids and meteorites provide key evidence on the formation of planetesimals in the 
Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation 
within a population of initially km-scale planetesimals. However, new models challenge this 
idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the 
gravitational collapse of small particles which have organised themselves in dense filaments 
and clusters in the turbulent gas. Particles concentrate passively between eddies down to the 
smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The 
streaming instability causes particles to take an active role in the concentration, by piling up in 
dense filaments whose friction on the gas reduces the radial drift compared to that of isolated 
particles. In this chapter we review new paradigms for asteroid formation and compare critically 
against the observed properties of asteroids as well as constraints from meteorites. Chondrules 
of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the 
primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled 
to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. 
We review recent progress on understanding the incorporation of chondrules into the asteroids, 
including layered accretion models where chondrules are accreted onto asteroids over millions 
of years. We highlight in the end ten unsolved questions in asteroid formation where we expect 
that progress will be made over the next decade. 


1. INTRODUCTION 

The Solar System contains large populations of pris¬ 
tine planetesimals that have remained relatively unchanged 
since their formation. Our proximity to the asteroid belt 
provides astronomers, planetary scientists and cosmo- 
chemists access to extremely detailed data about asteroid 
compositions, sizes and dynamics. Planetesimals are the 
building blocks of both terrestrial planets and the cores of 
the giant planets, as well as the super-Earths (with various 
degrees of gaseous envelopes) which are now known to or¬ 
bit around a high fraction of solar-type stars ^Fressin et al.\ 
2013) . The formation of planetesimals is thus a key step to¬ 
wards the assembly of planetary systems, but many aspects 
of the planetesimal formation process remain obscure. 

Recent progress in understanding planetesimal forma¬ 
tion was triggered by two important realisations. The first 
is that macroscopic dust particles (mm or larger) have poor 


sticking properties. Laboratory experiments and coagula¬ 
tion models show that it is difficult to form planetesimals 
by direct sticking of silicate particles, most importantly be¬ 
cause particle growth stalls at millimeter sizes where the 
particles bounce off each other rather than stick \Guttler 
20 lO). While some particle growth is still possible 


et al. 


at relatively high collision speeds, due to the net transfer of 
mass from small impactors onto large targets \Wunn efal. 


2005 Windmark et al. 2012a|l, the resulting growth rate i 


IS 


too low to compete with the radial drift of the particles. 

The second important realisation is that particles are con¬ 
centrated to very high densities in the turbulent gas flow. 
This idea is not new - Whipple ( 1972| l already proposed 
that large-scale pressure bumps can trap particles, as their 
radial drift speed vanishes in the pressure bump where the 
radial pressure gradient is zero. However, the advent of su¬ 
percomputing led to the discovery and exploration of a large 
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number of particle concentration mechanisms. Large-scale 

beyond asteroid formation (e.g. 

Cuzzi and Weidenschilling\ 

axisymmetric pressure bumps, akin to those envisioned by 

20061 Chiang and Youdin\ 2010 

\Johansen et al. 201411. 


in simulations of protoplanetary disk turbulence driven by 
the magnetorotational instability { Johansen et al.\ 2009a[ 
2012[). Particle densities reach high enough 


Simon et al. 


values inside these pressure bumps to trigger gravitational 
collapse to form planetesimals with sizes up to several 1000 
km Johansen et al.\ |2007| |201 ![ \Kato et aT\ |2012| l. The 
baroclinic instability, which operates in the absence of cou¬ 
pling between gas and magnetic field, leads to the forma¬ 
tion of slowly overturning large-scale vortices ^Klahr and\ 
Bodenheimer [2003 |l which can act as dust traps in a similar 


way as pressure bumps { Barge and Somme ria\ 1995| l. 

In the streaming instability scenario the particles play 
an active role in the concentration ^Youdin and Goodman\ 
2005^. The relative motion between gas and particles is sub¬ 


ject to a linear instability whereby axisymmetric filaments 
of a slightly increased particle density accelerate the gas 
towards the Keplerian speed and hence experience reduced 
radial drift. This leads to a run-away pile up of fast-drifting, 
isolated particles in these hlaments Johansen and Youdin 
|2007|l. The densities achieved can be as high as 10,000 


times the local gas density {Bai and Stone\pWiQ\ Johansen 


et al. 2012|l, leading to the formation of planetesimals with 


characteristic diameters of 100-200 km for particle column 
densities relevant for the asteroid belt, on a time-scale of 
just a few local orbital periods. 

A concern about large-scale particle concentration mod¬ 
els is that typically very large particles are needed for op¬ 
timal concentration (at least dm in size when the models 
are applied to the asteroid belt). Chondrules of typical 
sizes from 0.1 to 1 mm are ubiquitous in primitive mete¬ 
orites, but such small particles are very hard to concentrate 
in vortices and pressure bumps or through the streaming 
instability. One line of particle concentration models has 
nevertheless been successful in concentrating chondrules. 
Swiftly rotating low pressure vortex tubes expel particles 
with short friction times ^Squires and fsaton] |1990| |1991[ 


Wang and Maxey 1993|l. This was proposed to explain the 


The review is organised as follows. The first two sec¬ 
tions discuss the constraints on asteroid formation from the 
study of meteorites (Section]^ and asteroids (Section]^. In 
Section we review laboratory experiments and computer 
simulations of dust coagulation to illustrate the formidable 
barriers which exist to planetesimal formation by direct 
sticking. The turbulent concentration model, in which 
chondrule-sized particles are concentrated at the smallest 
scales of the turbulent gas flow, is discussed in Section 
The following Section|^is devoted to particle concentration 
in large-scale pressure bumps and through streaming insta¬ 
bilities. In Section]^ we discuss layered accretion models 
where the chondrules are accreted onto planetesimals over 
millions of years. Finally, in Section]^ we pose ten open 
questions in the formation of asteroids on which we expect 
major progress in the next decade. 

2. CONSTRAINTS FROM METEORITES 

Meteorites provide a direct view of the solid material 
from which the asteroids accumulated, while the crystalli¬ 
sation ages of the component particles and the degree of 
heating and differentiation of the parent bodies give impor¬ 
tant information about the time-scales for planetesimal for¬ 
mation in the Solar System. 

Meteorites may be broadly classified in two categories 
{Weisberg et al. 2006| l: primitive meteorites (also known 
as chondrites) and differentiated meteorites. Chondrites, 
which make up 85% of the observed falls, are basically 
collections of mm- and sub-mm-sized solids, little modi¬ 
fied since agglomeration and lithiflcation (compression) in 
their parent bodies. They exhibit nonvolatile element abun¬ 


characteristic sizes and narrow size ranges of chondrules 
observed in different chondrites \Cuzzi et ^|2001| l and lead 
to the formation of 100-km-scale asteroids from rare high- 
density concentrations ^Cuzzi ef Qr||2008||2010) l. However, 
the evaluation of the probability for such high-density con¬ 
centrations to occur over sufficiently large scales depends 
on scaling computer simulations to the very large separa¬ 
tions between the energy injection scale and the dissipation 
scale relevant for protoplanetary disks; Pan et al. ( 201 1[ ) 
found that the particle clustering gets less contribution from 
the addition of consecutively larger scales than originally 
thought in the model of \Cuzzi et a/.| ( |2008 2010| l. 

Therefore the incorporation of chondrules into the aster¬ 
oids is still an unsolved problem in asteroid formation. This 
is one of the main motivations for this review. We refer the 
readers to several other recent reviews on the formation of 
planetesimals which provide a broader scope of the topic 


dances comparable to the solar photosphere’s {Palme and 
Jones |2005l l. Differentiated meteorites derive from parent 
bodies which underwent significant chemical fractionations 
on the scale of the parent body, resulting in the asteroid¬ 
wide segregation of an iron core and silicate mantle and 
crust. In the process, differentiated meteorites have lost not 
only their accretionary (presumably chondritic) texture, but 
also their primitive chemical composition, for, depending 
on which part of the parent body they sample, some may be 
essentially pure metal (the iron meteorites) while others are 
essentially metal-free (the achondrites). 

It is among the components of chondrites that the oldest 


solids of the solar system, the refractory inclusions {Krot 

2005 |l, in particular Calcium- 


et al. 2004 MacPherson 


Aluminium-rich Inclusions (CAI), have been identified. 
Their age of 4567.3±0.16 Myr \Connelly et al. |2012| l 
marks the commonly accepted “time zero” of the Solar Sys¬ 
tem. But the ubiquitous components of chondrites are the 
eponymous chondrules {Hewins et al. 2005 [ Connolly and 
Desch\ 2004| l, which are millimeter-size silicate spherules 
presumably resulting from transient (and repeated) high- 
temperature episodes in the disk, but whose very nature 
remains a long-standing cosmochemical and astrophysical 
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2.1. Primary texture and aerodynamic sorting 

The texture of most chondrites has been reworked by 
impact fragmentation and erosion on their parent body. 
However, rare pieces of CM and CO chondrites have been 
found, referred to as “primary texture” ^etzler ef a/.[|199^ 


Brearley 19931, which seem to retain the nature of a pre- 


Onion Shell 


Rubble Pile 


Fig. 1.— The left panel indicates the initial state of a par¬ 
ent body of 100 km in radius, which was initially homoge¬ 
neous throughout, after heating by ^®A1. The proportions 
of the least (3) to most (6) altered material are constrained 
by meteorite statistics. The right panel indicates the same 
body after catastrophic fragmentation and reassembly, as a 
rubble pile with some highly altered material now near the 
surface. This body is then cratered by subsequent impacts, 
releasing samples of all metamorphic grades (adapted from 
\Scott and /?q/an||1981|l. 


brecciated body. Primary texture appears to consist of noth¬ 
ing but dust-rimmed chondrules of very similar properties, 
loosely pressed together. 

The constituents of most chondrites appear well-sorted 
by size, with strong mean size differences from one group 
to another ^Brearley and JonF^ 1998| l. Whether these dif¬ 
ferences arise from some regionally or temporally variable 
bouncing-barrier {Jacquet 2014a| l, some aerodynamic sort¬ 
ing process (Sections 5 and 7), or some aspect of the mys¬ 
terious chondrule formation process itself, they provide an 
important clue to primary accretion. Hezel et al. P008|l 


enigma [Boss 1996[ Desch et al. 20121. All these compo¬ 
nents are set in a fine-grained matrix comprised of amor¬ 
phous and crystalline grains native to the disk as well as 
rare presolar grains [Brearley 1996[ ). While composition- 
ally primitive, chondrites may have undergone some degree 
of thermal metamorphism, aqueous alteration and shock 
processing on their parent body. 

Despite their general petrographic similarity and roughly 
solar composition, chondrites are actually quite variable, 
and 14 distinct chemical groups have been recognized so 
far, each of which believed to represent a single parent body 
- sometimes with supporting evidence from cosmic-ray ex- 

200^ - or at least a family of 2.2. The abundance and distribution of ^®A1 and ®°Fe 

The melting of the parent bodies of differentiated mete¬ 
orites puts important constraints on the time-scale for plan- 
etesimal formation in the asteroid belt. While electromag¬ 
netic heating [Sonett and Colburn 1968| l or impact heat¬ 
ing [Kell et al. \991) have been considered in the litera¬ 
ture, the most likely source of planetesimal heating is the 
decay of the short-lived radionuclides (SLRs) ^^Al (with 
mean life-time t = 1.0 Myr) and ®°Fe (t = 3.7 Myr) (Urey 


have emphasized the need for better particle counting statis¬ 
tics, and indeed one recent chondrule size distribution mea¬ 
surement taken from Allende, of a far larger sample than 
analyzed previously ^Fisher et al.\ |2014| l, points to a distri¬ 
bution substantially broader for that chondrite than previ¬ 
ously reported. 

Aerodynamic sorting has been suggested often as an im¬ 
portant factor in selecting for the contents of primary tex¬ 
ture (see \ Cuzzi and Weidenschilling] |2006| for a review). 
Comparing the aerodynamical friction time of objects of 
greatly different density, such as silicate and metal grains, 
shows that their friction times are quite similar in the least 
altered meteorites [Dodd 1976[ Kuebler et al. 1999| l, sug¬ 
gesting that asteroids selectively incorporated components 
with specific aerodynamical properties (we discuss this fur¬ 
ther in Section|3.3|l. 


posure ages [Eugster et al. 


parent bodies formed in the same nebular reservoir (i.e. a 
compositionally distinctive space-time section of the disk). 
There are various levels of affinities between these groups 
(e.g. clans, classes) but we will be here content to distin¬ 
guish carbonaceous chondrites [henceforth “CCs”; com¬ 
prising the Cl, CM, CO, CV, CK, CR, CH, CB groups] from 
non-carbonaceous chondrites [henceforth “EORs”; com¬ 
prising the enstatite (EH, EL), the ordinary (H, L, EL) and 
the Rumuruti (R) chondrite groups]. Carbonaceous chon¬ 
drites are more primitive in the sense that they have a so¬ 
lar Mg/Si ratio and a ^®0-rich oxygen isotopic composition 
closer to that of the Sun (e.g. Scott and Krot\ 2003| l. Non- 
carbonaceous chondrites, though poorer in refractory ele¬ 
ments, are more depleted in volatile elements, have subsolar 
Mg/Si ratios and a more terrestrial isotopic composition for 


many elements (Trinquier et al. 20091. EORs have gener¬ 
ally undergone thermal metamorphism (see sketch in Eigure 
^ while aqueous alteration has been prevalent in CCs 


uss 


et al. 2006 1 Brearley 20031, but there are again exceptions. 


|1955[ ). Depending on their respective initial abundance, and 
on the time of planetesimal accretion, both could have sig¬ 
nificantly contributed to planetesimal heating. Additionally, 
short-lived nuclides provide crystallisation ages which can 
be calibrated using a long-lived radionuclide decay system 
such as Pb-Pb, under the assumption that the short-lived ra¬ 
dionuclide was homogeneously distributed in the solar pro¬ 
toplanetary disk. 

The content of SLRs in CAls is usually identified to 
that of the nascent Solar System ^auphas and Chaussi- 


don 


2011 1 . Excesses of ^®Mg linearly correlating with 


^^Al content were first observed in an Allende CAl in 1976 


[Lee et al. 1976|l. This isochron diagram demonstrated the 
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presence of in the nascent Solar System. The CAIs 
from a diversity of chondrite groups formed with an ini- 


tial (^®Al/^^Al)o of roughly 5x10^ (MacPherson et al. 
2 ^. 


A remarkably tight isochron for CAIs in the CV 

The 


( 20081 . 


chondrites was obtained by Jacobsen et al. 
deduced (^®Al/^^Al)o ratio of (5.23±0.13)x 10~^ is often 
considered as the initial value for the Solar System and the 
small dispersion as indicative of a narrow formation inter¬ 
val (< 40,000 yr). However, this interpretation should be 
limited to the region where the unusually large CV CAIs 

2009|l. This is especially true 


have formed (Krot et al. 


since many CAIs are known to have form ed without any 
^®A1 ( Liu et al. 2012 Makide et al. 2013| l. This indicates 
some level of heterogeneity in the ^®A1 distribution within 
the region where CAIs formed (assuming that region was 
unique, which is supported by the ubiquitous enrich¬ 


ment of CAIs compared to e.g. chondrules). Larsen et al. 
( |201 l| l used bulk magnesium isotopic measurements to sug¬ 
gest that the heterogeneity of ^®A1 distribution might have 
reached 80% of the canonical value in the solar protoplane¬ 
tary disk. However, Kita et al. (20131 wd \Wasserburg et'oT 
( 2012| l argue that the observed variations can be better as¬ 
cribed to small heterogeneities in the stable isotope ^®Mg. 

Although the presence of live ®°Fe in the early Solar Sys¬ 


tem was demonstrated almost 20 years ago (Shukolyukov 


and Lugmair 1993|l, the determination of the Solar Sys¬ 


tem initial abundance is complicated by the difficulty of 
obtaining good isochrons for CAIs \Quitte et al.\ |2007| l, 
given their low abundance in Ni. To bypass that difficulty 
most measurements were performed on chondrules which 
are believed to have formed from around the same time as 
CAIs up to three million years later ^Connelly et al. \ 2012| l. 
High initial (®°Fe/®®Fe)o ratios were originally reported in 
chondrules from unequilibrated ordinary chondrites (UOC) 
which experienced very little heating or aqueous alteration 


Solar System was lower than 1 x 10“®, though high levels 
of ®°Fe (up to 10“®) have been detected in some compo¬ 
nents of chondrites. At the time of writing this review there 
was no evidence of ®®Fe heterogeneity within the solar pro¬ 
toplanetary disk. 

2.3. The origin of ^®A1 and ®®Fe 

The initial Solar System ratio (^® AF^^Al)o of roughly 5 
X 10“® is well above the calculated average Galactic abun¬ 
dance ^uss et oTI |2009| l. This elevated abundance indi¬ 
cates a last minute origin for ^®A1. Production of ^®A1 by 
irradiation has been envisioned in different contexts (e.g., 
|Lee| |1978[ \Gounelle et al.\ |2006| l, but fails to produce 
enough ^®A1 relative to ^°Be (t = 2.0 Myr), another SLR 
whose origin by irradiation is strongly supported by exper¬ 
imental data \Gounelle et a7[[2^013| l. This leaves stellar de¬ 
livery as the only possibility for ^®A1 introduction in the 
Solar System. Though Asympto tic Giant Branch Star s pro¬ 
duce elevated amounts of ^®A1 (Lugaro et al. 2012 i, it is 
extremely unlikely that stars at this evolutionary stage are 
present in a star forming region (Kastner and Myers 1994| l. 
Thus massive stars are the best (unique) candidates for the 
origin of the Solar System’s Al. 

The recently obtai ned l ower estimates for the abundance 
of ®®Fe (see Section 2.2 1 are compatible with a Galactic 
background origin, independently of whether this is cal- 

p009] l, or 


culated crudely using a box model (Huss et al. 


taking into account the stochastic nature of star formation 
in molecular clouds \Gounelle et al. \ |2009| l. However, if 
the initial ®®Fe abundance corresponds instead to a much 
higher (®®Fe/®®Fe)o ratio of « 10“®, a last minute origin 
is needed. Irradiation processes cannot account for ®®Fe 


(e.g., Tachibana and Huss 2003 Mostefaoui et al. 2005 


production, because of its richness in neutrons (Lee et al. 
|I998| ). The winds of massive stars can also be excluded, 
given their low abundance in ®®Fe. In contrast, supernovae 
are copious producers of ®®Fe, essentially because this SLR 


is synthesized in abundance during the hydrostatic and ex¬ 
plosive phases of such massive stars \Woosley and HegTr 


20071. Two astrophysical settings have been envisioned so 


\Tachibana et al. 2006| l. Telus et al ( 2012| l showed that 
most of these previous data were statistically biased and 
that most chondrules do not show any ®®Ni excesses in¬ 
dicative of the decay of ®®Fe. The high values obtained 
in older publications could be due to statistical biases re¬ 
lated to low counts \Telus et aL\ |2012[ ) or to thermal meta¬ 
morphism which would have lead to the redistribution of 
Ni isotopes ^Chaussidon and BarrQf[ |2b09| l. Recently im¬ 
proved techniques for measuring bulk chondrules in un¬ 
equilibrated ordinary chondrites have yielded initial Solar 
System values for (®®Fe/®®Fe)o of «lxl0“® (Tang and 
\Daupha^\2Q\2\\Chen et ar||2013| l. This value is consis- disruptive to their environment (Chevalier 20001. 
tent with that inferred from Fe-Ni isotope measurements of 


far. In the first (classical) model ^Cameron and Truran\ 
|1977[ ), the supernova ejecta hits a molecular cloud core and 
provokes its gravitational collapse ^Boss and Reiser 2013| l 
as well as injecting ®®Fe and other SLRs. A newer model 
injects ®®Fe in an already formed protoplanetary disk (Hes- 
\ter et al.\ |2004t \Ouellette et al.\ |2007[ ). In either cases, the 
supernova progenitor mass is in the range 20-60 Mq, be¬ 
cause too massive supernovae are extremely rare and very 

Alter¬ 


a diversity of differentiated meteorites (Quitte et al. 2011 
Tang and Dauphas |2012|l. 


In conclusion, it seems that most CAIs formed with an 
initial ratio (^®AF^^Al)o of roughly 5 x 10“®, which can 
be considered in a first approximation as the Solar System 
average or typical initial value. Some heterogeneity was un¬ 
doubtedly present, but its exact level is still unknown. On 
the other hand, it is likely that the initial (®®Fe/®®Fe)o of the 


natively, generations of supernovae could have enriched the 
gas of the giant molecular cloud, so that the solar protoplan¬ 
etary disk simply inherited the elevated abundances of the 
birth cloud \Vasileiadis et q1.||2013| |. 

All these models suffer from an important difficulty, 
namely the overabundance of ®®Fe in supernova ejecta rela¬ 
tive to Al. The Solar System ®®Fe/^® Al mass ratio was ei¬ 
ther 4 X 10“® or 0.4 depending on the adopted initial ®®Fe 
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abundance. In any case this is well below the ®°Fe/^®Al 
mass ratio predicted by supernovae nucleosynthetic mod¬ 
els whose variation domain extends from 1.5 to 5.5 for a 
progenitor mass varying between 20 and 60 M©. Hetero¬ 
geneity in the composition of supernova ejecta has been 


proposed as a possible solution to that discrepancy {Pan 
\et al.\ |2012| l. However that variability is limited to a fac¬ 
tor of 4 which would help resolve the discrepancy only 
marginally in the case of the (unlikely) high ®°Fe value. Fi¬ 
nally, the astrophysical context of any supernova model is 
difficult to reconcile with observations of star-forming re¬ 
gions ^Williams and Gaidos\pnffT\\Gounelle and Meibom\ 
|2008| l. The commonly proposed setting for supernova con¬ 
tamination of a protoplanetary disk or or a dense core is 
similar to that of the Orion Nebula where disks are seen 
within a few tenths of parsec of the massive star 6^ C Ori. 
The problem with that setting is that, when 9^ C Ori will ex¬ 
plode in 4 Myr from now, these disks will have long evap¬ 
orated or formed planets. New disks or cores will have 
obviously formed by then but they will be at the outskirts 
of the 10 pc wide HII region created by 0^ C Ori ^Cheva- 


lier 


1999 1 . At such a distance, the quantity of ^®A1 deliv¬ 


ered into these disks or cores is orders of magnitude lower 
than the quantity present in the Solar System ^Looney et al. 
2006^. In other words, supernova remnants nearby dense 


phases are extremely rare {Tang and Chevalier 2014| l. In 
conclusion, it seems very unlikely that a nearby supernova 
was close enough to the Solar System to provide the known 
inventory of Al. 

If the low value of ®°Fe inferred from chondrites is cor¬ 
rect, then the presence of ^®A1 remains to be explained. 
Supernovae can be excluded as they would vastly overpro¬ 
duce ®°Fe (see above). The winds of mas sive stars have 
long been known to be Al and ^°Fe-poor ^rnould et al. 


1997|l. In the model of Gaidos et al. ( 2009| l and IToMwg 


( 2014| , ^®A1 is injected at the molecular cloud phase by 
the winds from a large number of Wolf Rayet stars. The 
problem with these models is that the Wolf Rayet phase 
is followed by the supernova explosion and therefore they 
produce a large excess of ®°Fe relative to Al and their re¬ 
spective Solar System values. To escape that caveat, [FoMwg 
( 2014| l has argued that Wolf Rayet stars do not explode as 
supernovae but directly collapse into black holes. Though 
this possibility has been theoretically envisioned, the recent 
observation of a Wolf Rayet star going supernova shows this 
is far from being the rule {Gal-Yam et al. 2014[ ). In addi¬ 
tion, models considering injection at the global molecular 
cloud phase cannot acco unt for the observed heterogeneity 
of Al (see Section 2.2 1 . 


The last class of models considered injection at the scale 
of single massive stars. Tatischejf et al. (20101 have en¬ 


visioned a single star that escaped from its parent cluster 
and interacted with a neighboring molecular cloud, inject¬ 
ing Al through its dense wind. It is compromised again by 
the proximity of the Wolf Rayet phase with the supernova 
phase. In addition. Wolf Ra yet stars are rare. In contrast 
Gounelle and Meynet (20121 proposed that ^®A1 has been 


injected in a dense shell of mass ~1000 Mq collected by 
the wind of a massive star. Evolutionary models of rotating 
stars are used, so the injection in the shell starts as early as 
the entry of the star onto the main sequence, lasts for some 
Myr and ends well before the supernova explosion. When 
the collected shell has become dense enough and gravita¬ 
tionally unstable, it collapses and a second generation of 
stars form which contain ^®A1. Detailed calculations have 
shown that as long as the parent star, baptized Coatlicue, is 
more massive than Mmi„ = 32 Mq, the abundance of ^®A1 
in the shell is equal to or larger than that of the Solar Sys¬ 
tem, depending on the mixing efficiency of the wind ma¬ 
terial with the shell. Because the mixing time-scale of the 
dense shell is comparable to its collapse time-scale, a cer¬ 
tain level of Al heterogeneity is expected ( Gounelle and 


Meynet 2012|l in agreement with observations. This model 


is in line with observations of induced star-formation within 
dense shells around massive stars {Deharveng et al. I IMoj i. 
Because it corresponds to a common - though not universal 
- mode of star formation, it implies that the Solar System 
is not the only of its kind to have formed with ^®A1, and 
that early differentiation of planetesimals might have been 
common in exo-planetary systems. |GoM»el/e] ( |2014[ ) has es¬ 
timated that the occurrence of planetary systems which are 
rich in ^®A1 and poor in ®°Fe is on the order of 1%. 

2.4. Timing of planetesimal accretion 

Planetesimal accretion itself cannot be dated directly 
with radioisotopic systems, since the mere agglomeration 
of different solids incurs no isotopic rehomogenization be¬ 
tween different mineral phases. Thus one can only obtain 
upper limits with the age-dating of pre-accretionary com¬ 
ponents (for chondrites) and lower limits from that of sec¬ 
ondary (parent body) processes. 

Important and ever-improving constraints have emerged 
since the publication of Asteroids III. In particular, Hf-W 
systematics of achondrites and irons have evidenced early 
differentiation, sometimes contemporaneous (within errors) 

[20T2| i. 


with refractory inclusion formation {Kruijer et al. 


This indicates that planetesimal formation started very early 


in the evolution of the solar system. Intriguingly, Libourel 


and Krot (2007 1 ascribed some olivine aggregates in chon¬ 


drites to this first generation of planetesimals (but see What- 
\tam et aL\ |2008[ \Jacquet et al.\ |2012a| for the alternative 
view that these formed in the protoplanetary disk). The 
above evidence for early differentiation is consistent with 
thermal modeling expectations, as the initial abundance of 

Al was sufficient to melt planetesimals, so that chondrites 
had to be accreted later to be preserved to the present day. 

Lower limits on chondrite accretion ages may be ob¬ 
tained from phases precipitated during aqueous alteration 


(see chapter by Krot et al). Fujiya et al. (20131 obtained 
ages of 4562.5-4563.8 Myr in Cl and CM chondrites (re¬ 
call the accepted age of CAI formation of 4567.3±0.16 
Myr, Connelly et aL \ 2012| l. As to non-carbonaceous chon¬ 
drites, recent Mn-Cr dating of fayalite formed during in- 
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cipient fluid-assisted metamorphism in the least metamor¬ 
phosed LL ordinary chondrites yields an age of 4564.9l|^'g 
Myr ^Doyle et al.\ 2014| l. Al-Mg systematics in the mildly 
metamorphosed H chondrite Sainte Marguerite indicate an 


age of 4563.1±Q.2 Myr {Zinner and Gdpel 20021. Ther¬ 
mal modeling based on ^^Al heating of the H chondrite 
parent body (one of the ordinary chondrite groups) con¬ 
strained by dating of chondrites of different metamorphic 
degrees indicate accretion ages ranging from 1.8 to 2.7 Myr 

and refer- 


2014 


after refractory inclusions (see Gail et al. 
ences therein), while [Ftpiya et al. P013 i advocate accre¬ 
tion of Cl and CM chondrites 3-4 Myr after “time zero” 
based on similar calculations. If ^®A1 really is the heat 
source behind chondrite alteration, then it would indeed 
make sense that carbonaceous chondrites, little affected 
by thermal metamorphism, accreted later than their non- 
carbonaceous counterparts [Grimm andMcSween 1993| l. 

Upper limits are provided by the ages of chondrules. In 
both carbonaceous and non-carbonaceous chondrites, chon- 
drule Pb-Pb ages in individual meteorites span a range of 
^4564-4567 Myr [Connelly et al. 2012| l, i.e. 0-3 Myr af¬ 
ter refractory inclusions, with younger ages reported by Al- 
Mg dating for CR ( Kita and Ushikubo] |2012| l and enstatite 
chondrite ^Guan et al 2006| l chondrules. Little correlation 
is seen between the age of chondrules and their composi- 
2012| l. The significance of Al-Mg 


tion ( Villeneuve et al. 


ages for chondrules, in particular in relationship with ap¬ 
parently somewhat older Pb-Pb ages (e.g. \Connelly et al.\ 
|2012^ , nevertheless remains uncertain. 

To summarize, it seems clear that differentiated planetes- 
imals accreted from the outset of the disk evolution while 
the known chondrites accreted later in the evolution of the 
disk (~2-4 Myr after refractory inclusions), as presumably 
required to escape differentiation, but the exact chronology 
of chondrite formation and alteration is still in the process 
of being firmly established. 

2.5. Accretion of chondrules directly after formation? 


Alexander et al. (20081 proposed that the retention of 
sodium, a volatile element, in chondrules during their for¬ 
mation indicated high solid densities in the chondrule- 
forming regions, up to 7 orders of magnitude above ex¬ 
pectations for the solar protoplanetary disk and possibly 
gravitationally bound (see also Cuzzi and Alexander 2006| l, 
although quite in excess of what the proportion of com¬ 
pound chondrules suggest ( Ciesla et al. 2004| l. This raises 
the possibility that the formation of chondrules and chon¬ 
drites, respectively, may have been contemporaneous, as 
also advocated by Metzler ( 2012| l based on the existence 
“cluster chondrites” comprised of mutually indented chon¬ 
drules (but see \Rubin and Brearley\\\99^ for a criticism of 
such hot accretion models). 


A further argument put forward by Alexander and Ebel 
( 2012| l is that chondrule populations in different chondrite 
groups are quite distinct ^Jones 2012| l. Indeed [Cuzz; et ai 
(|20I0|l noted that two populations of particles formed si¬ 


multaneously at 2 and 4 AU would be well-mixed within 1 
Myr for a = 10“^. Here a is the non-dimensional measure 
of the turbulent viscosity and diffusion. Relevant values 
for protoplanetary disks are discussed further in Section]^ 
but we note here that a value of 10“^ corresponds to the 
conditions which are expected if the mid-plane in the as¬ 
teroid formation region is stirred by turbulent surface layers 
(Oishi et al. 2007[ ). This would suggest that chondrules had 
to accrete rapidly to avoid homogenization. However, we 
do not know exactly the turbulence level or original space- 
time separations between chondrule-Zchondrite-forming lo¬ 
cations, as the asteroid belt may have undergone significant 
reshuffling ^Walsh ef^. 20111. Chondrite groups vary sig¬ 
nificantly in bulk composition. This indicates that there has 
been no thorough mixing of chondrite components, what¬ 
ever their individual transformations in the intervening time 
were, over the whole chondrite formation timescale. So 
whatever fraction of that time the period between chon¬ 
drule formation and chondrite accretion actually represents, 
it is not to be expected that chondrules should have been 
well-mixed over the different chondrite-forming reservoirs 
^acquet\\2.f)\A\)) . 

A link between chondrules and matrix is suggested, in 
the case of carbonaceous chondrites, by complementarity: 
the bulk meteorite is solar in some respect (e.g., the Mg/Si 
ratio) but its separate components (chondrules / matrix / 
refractory inclusions) are not (chondrules have a typically 
higher Mg/Si than solar while the converse is true for the 
matrix, \Hezel and Palme 2010| l. Complementarity - if 
verified, as for at least some elements it may reflect ana¬ 
lytical biases or parent body processes [Zanda et al. 2012} 
- would indicate a genetic relationship between chondrules 
and matrix, which would have exchanged chemical ele¬ 
ments upon formation, a relationship which would not be 
predicted e.g. in an X-wind scenario for chondrule forma¬ 
tion \Hezel and Palme\\20\0\ in which chondrules and ma¬ 
trix would have formed in different locations. But it does 
not require immediate accretion of chondrule and matrix. It 
only requires chondrules and dust grains to have been trans¬ 
ported in a statistically similar way, as was likely the case 
for a large portion of the disk evolution until accretion. Sev¬ 
eral batches of chondrule-ndust may have contributed to a 
given chondrite-forming reservoir provided again they suf¬ 
fered no loss of chondrules relative to dust or vice-versa 


[Cuzzi et al. 2005 Jacquet et al. 2012bI. This nonetheless 
does assume that at the stage of chondrule / matrix agglom¬ 
eration, there was no bias for or against the incorporation of 


any component [Jacquet 2014aI, which may be an impor¬ 


tant constraint on the accretion process. The problem is that 
small dust grains are much harder to incorporate into aster- 
oidal bodies than the macroscopic chondrules, due to their 
strong frictional coupling with the gas. One could never¬ 
theless envision that chondrules and matrix agglomerated 
together as compound objects (see Section 4.2.5) prior to 
incorporation in asteroidal bodies and/or that matrix-sized 
dust coagulated with ice into lumps with aerodynamical 
properties equivalent to chondrules. 
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Fig. 2.— The cumulative size distribution of asteroids N(>D), as a function of asteroid diameter D, from Morbidelli 


et al. p009|l. These coagulation models started with either km-sized planetesimals (left plot) or an initial size distribution 


following the current, observed size distribution of asteroids between 100 and 1000 km in diameter (right plot). The grey 
line shows the current size distribution of asteroids larger than 100 km in diameter. The model with small planetesimals 
overproduces asteroids smaller than 100 km in diameter (the upper dashed line represents the current size distribution of 
small asteroids while the lower dashed lines indicates a tighter constraint on the size distribution directly after accretion 
of the main belt). Starting with large asteroids gives a natural bump in the size distribution at 100 km in diameter, as the 
smaller asteroids are created in impacts between the larger primordial counterparts. 


So the jury is still out on whether chondrule formation 
immediately preceded incorporation in a chondrite or not. 
Given the chondrule age spread of 3 Myr within individual 
chondrites ( Connelly et al.\ 2012| l, as well as the presence 
of refractory inclusions and presolar grains which would 
not have survived chondrule-forming events, it is possible 
that chondrite components did spend up to a few Myr as 
free-floating particles in the gaseous disk prior to accretion. 


3. CONSTRAINTS FROM THE ASTEROID BELT 

The modern asteroid belt contains only a fraction of its 
original planetesimal population. However, the shape of the 
size distribution of the largest asteroids is primordial and 
gives important insights into the birth sizes of the planetesi¬ 
mals. Asteroid families provide a way to probe whether the 
asteroids are internally homogeneous or heterogeneous on 
large scales. 


3.1. Asteroid size distribution 


The observed size distribution of asteroids in the main 
belt shows a quite steep slope for bodies with diameter 
I?>100 km and a much shallower slope for smaller bod¬ 
ies (Bottke et al. 2005 | l. A similar change of slope with an 
elbow at i9~130 km is observed in the Kuiper belt popula¬ 
tion ^Fraser et a(.||2014l l. 

It was expected that the transition from a steep to a shal¬ 
lower slope is the consequence of the collisional disruption 
of smaller bodies. However, [.Bofffe et q(.| ( | 2005] | reached the 
opposite conclusion by examining the collisional evolution 
of the asteroid belt in detail. They used a number of con¬ 
straints (the total number of catastrophic asteroid families. 


the survival of the basaltic crust on Vesta, the existence of 
only 1-2 major basins on that body, etc.) to conclude that 
the integrated collisional activity of the asteroid belt had to 
be less than the one of the current main belt population in a 
putative time-span of 10 Gyr. If one supposes that initially 
the asteroid belt size distribution had a unique slope (the 
slope now observed for Z3>100 km), such a limited colli¬ 
sional evolution is not sufficient to reduce the slope of the 
size distribution of objects smaller than 100 km down to the 
observed value, i.e. to create the observed elbow. Therefore, 
\Bottke et a/.| ( |2005T ) concluded that the elbow at D~100 km 
is a fossil feature of the primordial size distribution. For 
the Kuiper belt, the constraints on the integrated collisional 
activity are not as tight as for the asteroid belt. Neverthe¬ 
less, models seem to suggest that collisional evolution alone 
could not create an elbow at diameters larger than 80 km 
{Pan and Sari 2005| l, which is significantly smaller than 
the observed value. 

Morbidelli et al. (2009 1 failed to produce the elbow at 
i9~100 km in the asteroid belt in collisional coagulation 
simulations starting from a population of small planetesi¬ 
mals (see Figure |^. So, having in mind the new models 
of formation of large planetesimals from self-gravitating 


clumps of chondrules or larger pebbles and boulders {Jo¬ 


hansen et al. 2007[ Cuzzi et al. 2008)1, they proposed that 


100 km was the minimal diameter of the original plan¬ 
etesimals. Moreover, not being able to reproduce the cur¬ 
rent slope of i9>100 km asteroids by mutual collisions be¬ 
tween bodies of 100 km in size, \Morbidelli et q/.|P 009| | 
argued that these large planetesimals were born with a sim¬ 
ilar slope. However, as we will see in Section]^ the cur¬ 
rent slope can be reproduced by considering the accretion 
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of chondrule-sized particles by 100-km-scale planetesimals 
during the gaseous disk phase, a process not considered by 
\Morbidelli et al. ( 2009| l. 


Weidenschilling (20111 managed to reproduce the elbow 


at Z3~100 km in the asteroid belt from collisional coagula¬ 
tion simulations starting from objects of 50 - 200 m radius. 
Because of the small size of these objects, collisional damp¬ 
ing and gas drag keep the disk very dynamically cold (i.e. 
with a small velocity dispersion among the objects). Hence, 
in the simulations Qf \Weidenschillin^ ( |201 l| l, the elbow at 
D~100 km is produced by a transition from dispersion- 
dominated runaway growth to a regime dominated by Ke- 
plerian shear, before the formation of large planetary em¬ 
bryos. However, any external dynamical stirring of the pop¬ 
ulation, for instance due to gas turbulence in the disk, would 
break this process. Moreover, these simulations are based 
on the assumption that any collision which does not lead to 
fragmentation results in a merger, but 100-m-scale objects 
have very weak gravity and the actual capability of bodies 
so small to remain bound to each other is questionable. Fi¬ 
nally we stress that the formation of 100-m-scale bodies is 
an open issue, in view of the bouncing barrier and meter- 
size barrier discussed in Section]?] 

3.2. Snowline problems 

Among the various meteorite types that we know, car¬ 
bonaceous meteorites (or at least some of them like Cl and 
CM) contain today a considerable amount (5-10%) of wa¬ 
ter by mass. Evidence for water alteration is wide-spread, 
and it is possible that the original ice content of these bod¬ 
ies was higher, close to the 50% value expected from un¬ 
fractionated solar abundances. Instead, ordinary chondrites 


contain < 2% water by weight (e.g./arosewic/t 1990 Krot 
et al.\ 2009 ]l; while ordinary chondrites do show signs of 


water alteration, it is unlikely that they ever contained 50% 
water ice by mass. This suggests that the parent bodies of 
these meteorites (Cl and CM vs. ordinary) formed on ei¬ 
ther side of the condensation line for water, also called the 
snowline. In other words, the snowline was located in the 
middle of the asteroid belt at the time when the asteroids 
formed. The Earth is also very water poor (some 0.1% by 
mass, although uncertain by a factor of ~5), suggesting that 
the planetesimals in its neighborhood were mostly dry and 
that water was delivered only by the small amount of plan¬ 
etesimals accreted from farther out ^Morbidelli et a/.] ] 2000 1 
Raymond et al. 2004] O’Brien et al. 2014) 1. The very dry 


enstatite chondrites have been proposed to arise from an ex¬ 
tinct portion of the asteroid belt located between 1.6 and 2.1 
AU from the Sun ^Bottke et aZ.]]2012l l. 

The problem with this picture is that chondritic bodies 
accreted late (2-4 Myr after CAIs) and that the snowline is 
expected to migrate towards the Sun with time. The temper¬ 
ature of the disk is set by the equilibrium between heating 
and cooling. In the inner part of the disk, the heating is 
predominantly due to the viscous friction of the gas in dif¬ 
ferential rotation around the Sun, so it is related to the gas 


accretion rate of the star. The cooling rate is governed by 
how much mass is in the disk in the form of micron-sized 
grains. In Bitsch et al. ( 20I4al l, the snowline is at 2-3 AU 
only in the earliest phases of disk evolution when the accre¬ 
tion rate is M = Mq yr“^; but the snowline migrates 
to 1-2 AU already after 500,000 yr when the accretion rate 
drops below M = 10~^Mq yr“^. These models are proba¬ 
bly optimistically warm, because it is assumed that the ratio 
between the mass in /rm-sized dust and the mass in gas is 
1% (this reflects the solar metallicity but planetesimal ac¬ 
cretion should eventually decrease the dust content) and ac¬ 
count for the heating produced by stellar irradiation by a 
star more luminous than our early Sun. According to the 
observation of the accretion rate of stars as a function of 
age ( Hartmann et al. 1998)1 and photoevaporation models 
Alexander and Armitage\ 2006] ), the accretion rate on the 
star should drop to M = IO^^Mq yr“^ at 2-3 Myr, i.e. at 
the time of chondrite formation. The fact that there is no 
sign of primary accretion in the Solar System after ~3 Myr 
also suggests that the disk disappeared at that time (or that 
all the solids could have been incorporated into parent bod¬ 
ies at that time - however, accretion of solids is unlikely to 
reach 100% efficiency). Thus, at the time of chondrite for¬ 
mation, the snowline should have been well inside the inner 
edge of the asteroid belt, possibly even inside 1 AU. This is 
at odds with the paucity of water in the ordinary chondrites 
and in the Earth. 

A way to keep the disk warm has been proposed recently 


by Martin and Livio (2012 20131. In their model, a dead 
zone in the disk, with low turbulent viscosity, becomes very 
massive and develops gravitational instabilities which heat 
the gas. The unlimited pile-up of gas in the dead zone, how¬ 
ever, could be a consequence of a simplified disk transport 
model, and is not observed in more complex 3D hydrody- 
namical calculations \Bitsch et a(.]]2014b| l. Also, the model 
of Martin and Livio { 2012) predicts a cold region where ice 
could condense inside the orbit of the Earth, which seems 
inconsistent with the compositions of Venus and Mercury. 

A possible way to keep the asteroid belt ice-poor, de¬ 
spite freezing conditions, is to stop the inwards flow of 
icy pebbles in a pressure bump. The inwards drift of peb¬ 
bles is caused by friction with the slower moving gas (this 
effect is discussed more in Section 16.1) l. Then when the 
nebula cooled locally below the ice condensation temper¬ 
ature, there was no water vapor present to condense, and 
the planetesimals formed there would have a volatile deple¬ 
tion similar to that achieved in a warm disk. \Bitsch et al.\ 
(2014b]) explored the effects of viscosity transitions at the 
snowline (as suggested by Kretke and Lin 2001) , but found 
that only very steep radial gradients in the a parameter al¬ 
low the formation of a pressure bump. \Lambrechts et al.\ 
(2014 1 showed that the formation of a proto-Jupiter of about 
20 Me can produce a strong pressure bump just exterior to 
its orbit. If proto-Jupiter formed at the snowline when the 
snowline was at the outer edge of the asteroid belt, then its 
presence would have shielded the asteroid belt and the ter¬ 
restrial planet region from the flow of icy pebbles, while 
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small silicate dust and chondrules would remain for longer 
times in the asteroid belt due to their slow radial drift. 

3.3. Are asteroids internally homogeneous? 

The possibility that the asteroids grew incrementally by 
layered accretion of chondrules (Section]^ implies that as¬ 
teroids are internally heterogeneous, in the most extreme 
case with a chondritic surface layer residing on a differenti¬ 
ated interior. We discuss here briefly whether such internal 
heterogeneity is supported by observations of the asteroid 
belt. \Burbine et a/.| ( [2002| note that 100-150 distinct me¬ 
teorite parent bodies, 3/4 of them differentiated, are repre¬ 
sented in the meteorite collection. However, this sample 
is biased towards the sturdy irons and against the weaker, 
never-melted primitive chondrites. 

Differentiation is a constraint on formation age. Most 
studies have suggested that sizeable asteroids forming < 
10® years after CAIs are almost certain to have thoroughly 
melted, but those which formed more than 2 Myr after CAIs 
may have escaped melting except near their centers. More 
recent work paints a more complicated picture (see chap¬ 
ter by Elkins-Tanton et ai). Specifically, the Allende parent 
body, source of primitive CV3 chondrites, is thought to have 
melted near its center, as evidenced by the paleomagnetism 
detected in these meteorites ^Elkins-Tanton et oT] 12011)1; 


other CV chondrites may show a similar signal (We/ii ef a/. a rubble pile (Figure 1, Taylor et al. 


sual and near-IR wavelengths ^Mothe-Diniz and Nesvorn^ 
|2008| ). This conclusion is nevertheless complicated by the 
identification of interlopers in (and exclusion from) the as¬ 
teroid families. Hence there is an inherent tendency to ob¬ 
serve low albedo variations within asteroid families. In con¬ 
trast, the Eos family and the Eunomia family have unusually 
large internal variance \Mothe-Diniz et a/[||2005[[2008[ ) and 
the Eunomia family looks like what an internally differen¬ 
tiated S-type might produce (see discussion in \Weiss arid 
\Elkins- Tanton\\2Q\3\ . 

The three ordinary chondrite groups (H, L, and LL) are 
each thought to derive from a single parent body, based on 
a clustering of cosmic ray exposure ages of chondrites in 
each group, as if they were excavated by the same few large 
impacts { Eugster et al. \ 2006b| l. The thermal history of the 
H chondrite parent body under internal heating by Al was 
modeled by Trielojf et al. ( 2003| l, Monnereau et al. ( 2013| l, 
and Henke et al. ( 2013| l who all concluded it was a roughly 
100-km-radius body. Chondrites from all three OC groups 
show different amounts of thermal alteration from this pro¬ 
cess, designated as metamorphic grades 3-6. Models im¬ 
ply that the H chondrite parent, at least, first incurred ther¬ 
mal alteration in an “onion-shell” fashion, with the most 
strongly heated H5/H6 material heated deep in the interior, 
and then was catastrophically disrupted and reassembled as 

1987) \Scott et 'ai. 


2010) . This suggests that the CV parent body was differen- |2014| ). The small fractional abundance of minimally al- 


tiated in its interior, but preserved an undifferentiated chon¬ 
dritic layer. Other authors have ascribed the magnetic field 
in the CV chondrites to impacts (e.g. Wasson et al. 2013| l. 

Modeling of the buoyancy of silicate melts of different 
composition suggests that, for C-type composition, melt 
might be dense and remain stable at depth, but for S-type 
(OC) compositions, and certainly for Enstatite composi¬ 
tions, melt is less dense than surrounding material and 
will rise to manifest on the surface ( |Fm and Elkins-Tanton\ 
2014). Thus there might be old, centrally-melted C-type as- 


tered chondrites (H3/H4; Eigure 1) constrains the accretion 
of the H parent body t o have happened quickly - probably 
faster than 3 x 10® yr {Henke et al.\ 2013 Vernazza et al. 


20141). The constant mass growth rate adopted in Henke 


et al. 


(20131 is nevertheless not applicable to layered accre- 


teroids such as the Allende parent body, where the evidence 
for differentiation remains buried, but the lack of evidence 
for significant surface melting on most S-type asteroids may 
argue that most of them remain unmelted and undifferenti¬ 
ated throughout (see chapter by Elkins-Tanton et al., and 
IVTem and Elkins-Tanton\\20l3) . 

We can sample the internal properties of asteroids in two 
ways - from the observed color and albedo distributions of 
collisionally disrupted asteroids in families (see also chap¬ 
ters by Nesvorny et al. and Michel et al.) and from the mete¬ 
orites that derive from them (see next paragraph). The NE¬ 
OWISE mission has measured the albedos of a large num¬ 
ber of asteroids and families, finding a dichotomy in albedo, 
roughly corresponding to the classical S and C types, that 
is most evident in the outer main belt. Collisional families 
cover the entire belt, so they avoid the sampling bias that af¬ 
fects meteorites. Many families have internal albedo distri¬ 
butions that are narrower than the global spread of albedos 
across families. A similar story is told by observations at vi¬ 


tion, which results in run-away accretion and deposition of 
most of the mass towards the end of the growth phase (see 
further discussion in Section]?). 

While the major element chemical compositions, and the 
oxygen isotopic compositions, of the OC groups differ sig¬ 
nificantly, there is little or no discernible variation of ei¬ 
ther chemical or isotopic composition across metamorphic 
grades in any of the three groups { Jarosewi^ 1990 \Wood\ 
20051; see also Tables 1 and 2 of \Clayton et al. (19911. The 


variation of chondrule size across metamorphic grade is less 
well studied (A. Rubin, personal communication 2014). 

^odd\ ( |1976) demonstrated that the difference in the 
aerodynamical friction time between metal grains and sil¬ 
icate chondrules can explain silicate-metal fractionation in 
the ordinary chondrites. In this picture the LL chondrites 
are under-abundant in metal because the parent body was 
successfully able to accrete large silicate-rich chondrules 
(the chondrules in the LL chondrites are larger than those 
in both H and L chondrites, see e.g. Nelson and Rubin\ 


20021. The oxidation of metallic Fe could have happened as 


the LL chondrite parent body accreted water-bearing phyl- 
losilicates \Rubin 2005) . In the layered accretion model 
presented in Section |7) asteroids will accrete larger and 
larger particles as they grow. Metal grains had a more re¬ 
stricted size range and hence a relatively small parent body 
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of the H chondrites would accrete mainly small chondrules 
together with metal grains. A small size of the H chon- 
drite parent body is further supported by the low fraction 
of strongly metamorphosed samples that we have from that 
group (Dodd 1976|l. A similar story is told for the enstatite 


chondrites [Schneider et al. 1998|l: the EH group has more 


metal and smaller chondrules than the EL group. This aero¬ 
dynamical size sorting may be evidence of asteroid growth 
by layered accretion (Section l7 1 or asteroid formation by 
turbulent concentration (Section 

4. DUST GROWTH BY STICKING 

Now that we have given an overview of some of the con¬ 
straints from meteorites and asteroids, we can turn to the 
theoretical models of planetesimal formation. In this sec¬ 
tion we discuss the growth of dust by direct sticking; sub¬ 
sequent chapters discuss gravitational instability models. 

4.1. Particle-gas interaction 

The dynamical behavior of a particle in gas depends on 
both its size and density, as determined by its friction time Tf 
in the nebula gas ( Whipple] \912\^eidenschilling 1977al l. 
Eor particles smaller than the gas molecule mean free path 
(approximately 10-100 cm in the asteroid belt region, de¬ 
pending on the uncertain value of the gas density), the fric¬ 
tion time is 


Tf = 


ap. 


( 1 ) 


where a and p, are the particle radius and density, and uth 
and pg are the gas thermal speed and mid-plane density. 
The thermal speed of the gas molecules is in turn connected 
to the isothermal sound speed through Uth = The 

Stokes number St is defined as St = f7rf, where 17 is the 
(local) orbital frequency of the protoplanetary disk. The 
translation from Stokes number to particle size follows 


(2/7r)rgSt 


80cmSt/g(r)(^^) ' ' . (2) 


Here ii'g = x/^pgiTg is the gas column density and fg{r) 
a parameter which sets gas depletion relative to the MMSN 
as a function of semi-major axis r. We set in the second 
equality p, = 3.5g/cm^, relevant for chondrules. The 
Stokes number controls many aspects of dust dynamics. 
Particles of larger Stokes numbers couple increasingly to 
larger, longer-lived, and higher-velocity eddies in nebula 
turbulence, thus acquiring larger relative velocities. Solu¬ 


tions for particle velocities have been developed by Voelk 


et al. 


( |1980| l and Mizuno et al. (1988|l, including closed- 


form analytical expressions by Cuzzi and Hogan\\2QQ3) and 
Ormel and Cuzzi (20071. Importantly, the Stokes number 


also controls the degree of sedimentation, with the scale- 
height of the sedimented mid-plane layer, i7p, given by 




( 3 ) 


Here Hg is the gas scale-height and a is a measure of the 
turbulent diffusion coefficient D normalised as D = aCsHg- 


{s&& Johansen et al. 2014 for references). Significant sedi¬ 


mentation can occur when St > a. Chondrules of mm sizes 
have a Stokes number of St ~ 10“^; hence chondrules will 
only settle out of the gas if a 10“^. 

4.2. Dust growth 

The study of dust growth has been an extremely active 
field, both experimentally and numerically, since Asteroids 


III. Subsequent reviews were presented by Dominik et al. 
(2007]l, \Blum and Wurrri\ (| 2008| l and Chiang and Youdin 
( 2010| l; the older review by Beckwith et al. ( 2000| l is also 
valuable for basics. Two very recent overview chapters pre¬ 
senting both the basic physics and selected recent highlights 
are by Johansen et al. ( 2014[ ) and Testi et al. ( 2014| l; for 
efficiency we build on those chapters and here emphasize 
specifics relevant to asteroid formation. 

4.2.1. Sticking 

Particles can stick if their relative kinetic energy ex¬ 
ceeds certain functions of the surface energy of the mate¬ 
rial, which depends on composition. At the low relative 
velocities for small monomers (0.1-10 pm) under nebula 
conditions, both ice and silicate particles stick easily and 
form loose, porous aggregates. The process continues until 
the aggregates are at least 100 pm in radius. Aggregates 
can continue to grow and stick at larger velocities, if their 


open structure is able to deform and dissipate energy (Wada 
2009|l. The entire process of growth beyond roughly 


et al. 


100 pm fluffy aggregates depends on just how much these 
aggregates can be compacted by their mutual collisions. 
Recent studies that concentrate on icy particles outside the 
snowline have argued that the high surface energy of ice 
prevents significant compaction from occurring (and keeps 
relative velocities small) until particles have grown to ex¬ 
tremely large size - hundreds of meters - with extremely 
low density \Okuzumi et al. 2012|l. 


4.2.2. Bouncing 

The surface energies of silicates are much smaller than 
those of ice, so it is easy for even mm-sized silicate parti¬ 
cles to compact each other in mutual collisions. Relative 
velocities large enough to cause compaction and bouncing 
are acquired by roughly mm-size silicate particles in nebula 
turbulence. Coagulation modeling by Zsom et al. |(|20T0ll, 
consistent with experiments ^Guttler et al. 2010[ Weidling 
\et al.\ |2012| l, revealed a bouncing barrier in this size range 
where growth of silicate aggregates by sticking ceased. This 
new banier joins the long-known fragmentation barrier and 
radial drift bander which, even if the bouncing barrier can 
be breached, tend to frustrate growth in the asteroid belt re¬ 
gion beyond dm-m size ^Brauer et al.\\2002i\^irnstiel et al.\ 

[20To1[2^ . 
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4.2.3. Fragmentation 

A simple critical velocity can be used to refer to 
fragmentation of two comparable masses. This approach 
has been modified in some treatments to include some mass 
transfer by a smaller projectile hitting a larger target at high 
velocity, even if the projectile is destroyed and some mass 
is ejected from the target ^Wurm et al.\ 2005) 1. An alternate 
approach is to treat the fragmentation threshold as a critical 
kinetic energy per unit mass Q*, which can be thought of 


as a critical velocity squared for same-size particles [Stew- 
\art and Leinhardi\ |2009| l. The latter treatment automati¬ 
cally accounts for particle size differences and thus allows 
accretion of small particles to proceed at collisio n veloci- 
ties much higher than the nominal Ufrag V^- Stewart 


and Leinhardt\{ 2009 1 treated solids as weak rubble piles, all 


calibrated using experimental work hy \Setoh et q1.| ( |2007| i. 
These expressions allow for the higher efficiency of low- 
velocity collisions in fragmentation than for hypervelocity 
impacts. For particles made of small silicate grains, a value 
of Q* on the order of 10"^ cm^/s^ is suggested, with an asso¬ 
ciated Ufrag ~1 to several m/s, for weak cm-m-sized aggre¬ 


gates (5'c/zrd/7/er ef a/. 2012 Stewart and Leinhurdt 2009 
\Wada et ~ 2009 2013 |l or even less ( Beitz et al. 20 


4.2.4. Lucky particles 

The bouncing barrier, in preserving most of the avail¬ 
able solids at small sizes, may provide a target-rich envi¬ 
ronment for growth of much larger “lucky” particles, which 
experience few, or low-velocity, collisions and avoid de¬ 
struction while steadily growing from much smaller parti¬ 
cles to large sizes ( Windmark et al. 2012a|b Garaud et al. 


201^ Dr^zkowska et al. 2013|l. Ho wever, Windmark et al. 
( 2012a|b| l and Garaud et al. ' p0T3] l did not include radial 
drift, which is important because the growth time for these 
lucky particles greatly exceeds the drift time to the Sun. 
Drqzkowska et al. ( |2013| l removed the radial drift prob¬ 
lem with a pressure bump at the edge of a dead zone, and 
still found that the total number of meter-size particles they 
could produce in 3 x 10^ years was in the single digits. An 
isolated particle cannot trigger collective effects such as the 
streaming instability (Section 6), but can only keep growing 
by sweep up. Johansen et al. ( 2008| l and Xie et al. ( 2010| ) 
have modeled this “snowball” stage and find that growth in 
this fashion is extremely slow unless the nebula turbulent a 
is very low, because the small feedstock particles are verti¬ 
cally diffused to a low spatial density otherwise. 


small aggregates {MacPherson et al. 1985 Metzler et al. 
\\992\\Blandetal.\\2Qn) . Matrix and rims were reviewed in 


depth by Scott et al. ( 1984| l; amongst many interesting re¬ 
sults, they found evidence that rims were accreted as numer¬ 
ous aggregates of variable mean composition, rather than 
as monomers. Rubin ( 201 1| ) suggested that carbonaceous 
chondrites formed in dusty regions of the solar protoplan¬ 
etary disk and that matrix accumulated into mm-cm-sized 
highly porous dust balls. In this picture, chondrules ac¬ 
quired matrix rims by collisions with these dust balls rather 
than in collisions with smaller particles. 

The accretion of fine grained rims was modeled by|Mo 


lor-\ 


fill et al. ( 1998 1, Cuzzi ( |2004 1 , and Ormel et al. ( |2008 1. In 
these models, relative velocities predicted for the nebula en¬ 
vironment have been shown to be compatible with sticking 
and compaction under the theory oi^ominik and Tielens 
( |1997| l; [OrOTe/ et al] ( |2008| l added the effects of interchon- 
drule collisions on compacting the rims further. Ormel et al. 
( 2008) 1 found that chondrules which had accreted porous 
rims of dust between collisions could stick more easily be¬ 
cause the porous rims acted as shock absorbers, resulting in 
composite, dm-size objects formed of rimmed chondrules. 

4.2.6. Gravitational scattering barrier 

There is one further barrier for particles trying to grow 
incrementally, by sticking, to km and larger sizes. This 
gravitational scattering barrier arises because, in turbu¬ 
lence, nebula gas has small density fluctuations associated 
with pressure and vorticity. These density fluctuations, pri¬ 
marily on large scales, scatter growing planetesimals to 


achieve high relative speeds \Laughlin et al. 

12004) Nel- 

son and Papaloizou\ 2004) Nelson and Gressel 

)2010) Yang 

et al. 2012)1. The random velocities acquired by 1-10-km- 


scale planetesimals in this way are sufficient to put them 
into an erosive, rather than accretionary, regime throughout 
most of the Solar System for a range of the most plausible 
a values [ Ida et al. )2008 \Stewart and Leinhardt\ )2009) l. 
Ormel and Okuzumi ( 2013) l found that planetesimals need 
to have radii of 100 km or more to be able to survive these 
destructive encounters. 

Overall, the barriers to formation of 100-km-scale aster¬ 
oids by incremental growth-by-sticking appear formidable. 
For these reasons, a number of models have arisen which 
avoid these problems with a leapfrog process in which 10- 
1000-km-scale asteroids form directly from smaller parti¬ 
cles stuck at one of these barriers. 


4.2.5. Chondrule rims and chondrule aggregates 

The fine-grained component of chondrites is not only 
found in a featureless background matrix; it is also found 
rimming individual chondrules and other coarse particles, 
often Ailing cavities [Metzler et al. 1992) Brearley 1993) 
Cuzzi et al. 2005)1. The origin of these^ne grained rims has 


been debated [Lauretta et al. 2006) 1. One school of thought 
regards them as accretionary rims, swept up as a cooled 
chondrule moves relative to the gas and entrained dust or 


5. TURBULENT CONCENTRATION 

The leapfrog model sometimes referred to as turbulent 
concentration (TC) or turbulent clustering can in princi¬ 
ple make 10-100-km-radius planetesimals directly from 
chondrule-sized particles. The turbulent concentration 
model was motivated originally by laboratory and numeri¬ 
cal experiments that showed dense particle clusters forming 
spontaneously in isotropic turbulence, in which the most 
intense clustering was seen for particles with friction time 
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Tf equal to the Kolmogorov (or smallest eddy) timescale 
{Squires^md^atot^ 1990 [1991 Wa«g and Maxey 1993 


\Eaton and Fessler 1994 1. It was realized that this condition 
was very closely satisfied under canonical nebula condi¬ 
tions by the very chondrule-size particles that make up the 


bulk of primitive chondrites ( Cuzzi et al. 1996 20011. 


5.1. A brief history of TC models 

What remains uncertain and a subject of current study, 
as with many properties of turbulence, is just how the labo¬ 
ratory experiments and numerical simulations translate into 
actual nebula conditions. The flow regime is described by 
its turbulence Reynolds number Re = UL/v = aCsHg/v 
where U and L are some typical velocity and length scale 
of the flow, and v is the molecular viscosity. All current 
experiments and simulations are run at Re far smaller than 

(2001)1 noticed that the 


likely nebula values. Cuzzi et al. 


Hogan and Cuzzi (20071 showed that mass-loading feed¬ 


back of the particle burden on gas turbulence caused the 
concentration process to saturate at a mass loading $ = 
Pp/Pg ~ 100’ precluding the small, dense clumps advo¬ 
cated by Cuzzi et al. P001| l. Based on this, Cuzzi et al. 
( 2008| l advocated gravitational binding and ultimate sedi¬ 
mentation of much larger clumps, on the order of 10^ — 10'* 
km, with 4) 10 — 100. They showed that a long-neglected 


number of caveats regarding the built-in assumptions of this 
model; one caveat regarding scale-dependence of the con¬ 


centration process has been found to be important enough 
to change the predictions of the scenario quantitatively (see 


below). Subsequently, Cuzzi and Hogan (20121 resolved a 


discrepancy in a key timescale between Cuzzi et al. ( 20T0| 
and Chambers ( 2010| l, which makes planetesimal formation 
1000 times faster than in Cuzzi et al. ( 2010| l (and corre¬ 
spondingly slower than m\Chambers 2010|l. 


volume fraction of dense clumps increased with increasing 
Re, and suggested a scaling which would map the behavior 
to the much higher values of Re relevant for turbulence in 
protoplanetary disks. 


5.2. New insights into turbulent concentration 

The primary issues are whether it is always Kolmogorov 
friction time particles that are most effectively concen¬ 
trated, and whether the physics of their concentration is 
scale-invariant. \Hogan and Cuzzi] ( |2007| ) argued by anal¬ 
ogy with the observed scale-invariance of turbulent dis¬ 
sipation, which is dominated by Kolmogorov-time vortex 
tubes (little tornados in turbulence) that the concentration 
of Kolmogorov-friction-time particles would also be scale- 
invariant (see also Cuzzi and //ogaw) ( 2012| l. They devel¬ 
oped a so-called cascade model by which to extend the low- 
Re results to nebula conditions. 


The primary accretion scenarios of Cuzzi et al. (2010 1 


and Chambers (2010 1 used this cascade model to gener¬ 
ate density-vorticity PDFs as a function of nebula scale. 
Pan et al. (|20]4Jl ran simulations at higher Re than |77oga»| 


and Cuzzi (2007 1 and found that the clump density PDFs 


dropped faster than would be predicted by the scale- 
invariant cascade. They suggested that the physics of parti¬ 
cle concentration might indeed be scale-dependent, and that 


discovery by Sekiya (19831 precludes genuine, dynamical- 
timescale collapse for dense clumps of chondrule-size parti¬ 
cles under plausible conditions. Particles of these sizes and 
friction times can only sediment slowly towards the center 
of their bound clump, on a much longer timescale (10^ — 10^ 
orbits). \Cuzzi et ar| ( |2008| l suggested a criterion for stabil¬ 
ity of these bound, yet only slowly-shrinking, clumps which 
led directly to the conclusion that they would preferentially 
form large planetesimals comparable in size with the mass- 
dominant 100 km radius mode advocated for the early as¬ 
teroid belt by \Bottke et al.\ ( |2005[ ) - thus leaping over the 
long-troublesome mm-m-km-size barriers. Objects formed 
in this way are compatible with the primary texture seen 
in the most primitive unbrecciated CM and CO chondrites 
^Metzler et al. 1992) ^rearley\ 1993) 1, and would display 
deep homogeneity of bulk chemical and isotopic properties. 

\Cuzzi et ar\\20\0\ and )C/iamher5) ( )20T0) ) went on to de¬ 
scribe an end-to-end primary accretion scenario, combin¬ 
ing stability thresholds with calculated probability distri¬ 
butions of clump density, finding that a range of nebula 
conditions (all implying > lOx local enhancement of the 
usually assumed 1% cosmic solids-to-gas ratio within some 
few 10*^ km of the midplane) could match the required rate 
of planetesimal formation and the characteristic mass mode 
around 100-200 km diameter. Cuzzi et al. (20101 gave a 


planetesimal formation rates obtained using the Hogan and 
\Cuzzi\^iXfl\ cascade might be significantly overestimated. 

Ongoing work supports this concern about scale depen¬ 
dence. Cuzzi et al. ( 2014)l hav e analyzed much higher Re 
simulations [ Bee et al. 2010) 1 and found that the cascade 
measures called “multiplier distributions”, which determine 
how strongly particles get clustered at each spatial scale, do 
depend on scale at least over the largest decade or so of 
length scale; that is, the scale-invariant inertial range for 
particle concentration and dissipation does not become es¬ 
tablished at the largest scale, causing little concentration to 
occur until roughly an order of magnitude smaller scale. 
Because the cascade process is multiplicative, this slow start 
means that fewer dense zones are to be found at any given 
scale size than previously thought. 

New, scale-dependent cascades can now be implemented 


using the results of Cuzzi et al. (2014 1 at the largest scales. 
The quantitative implications are not clear as of this writing. 
One qualitative result is that particles with friction times 
2-10 X longer than originally suggested are more strongly 
clustered at the most relevant (larger) length scales (con¬ 
sistent with previous work by Bee et ai\ 2007) Pan et al. 
)201 1) ). Because we know the size and density of chon- 
drules, the implications will be different preferred values 
of gas density and/or turbulent a (see Cuzzi et al. )20mT l. 
A related implication is that concentrated particle sizes 


may not be as narrowly peaked as shown in Cuzzi et al. 


()2001)l. However, at the same time, new observations are 
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also showing broader chondrule size distributions in chon¬ 
drites fisher et a/.||2014| l. 

6. PRESSURE BUMPS AND STREAMING INSTA¬ 
BILITY 

The turbulent concentration mechanism described in the 
previous section operates on the smallest scales of the turbu¬ 
lent flow (although the vortical structures which expel par¬ 
ticles can be very elongated). The dynamical time-scales 
on such small length-scales are much shorter than the local 
orbital time-scale of the protoplanetary disk. In contrast, 
the largest scales of the turbulent flow are dominated by the 
Coriolis force, and this allows for the emergence of large- 
scale geostrophic structures (high-pressure regions in per¬ 
fect balance between the outwards-directed pressure gradi¬ 
ent force and the inwards-directed Coriolis force). 


Whipple (19721 found that particles are trapped by the 


zonal flow surrounding large-scale pressure bumps. Pres¬ 
sure bumps (in a way azimuthally extended analogs to the 
vortices envisioned in Barge and Sommen^ 19951 can arise 
through an inverse cascade of magnetic energy {Johansen 


\et al.||2()09at|5»wow et al.\\2Q\2\^ittrich et ar||2013| l in tur- 

bulence driven by the magnetorotational instability ^Balbus 
\and Hawley\\\99\\ . Pressure bumps concentrate primarily 
large (0.1-10 m) particles which couple to the gas on an or¬ 
bital time-scale ^Johansen et al.\ 2006| l, reaching densities 
at least 100 times the gas density which leads to the forma¬ 
tion of 1000-km-scale planetesimals {Johansen et al. 2007 


|2011| l. The magnetorotational instability is nevertheless no 
longer favoured as the main driver of angular momentum 
transport in the asteroid formation region of the solar pro¬ 
toplanetary disk, since the ionisation degree is believed to 
be too low for coupling the gas to the magnetic field (see 
review hy\Turner et al. \ 2014|l. 


The magnetorotational instability can still drive turbu¬ 
lence (with a in the interval from 10“^ to 10“^) in the mid¬ 
plane close to the star (within approximately 1 AU where 
the ionisation is thermal) and far away from the star (beyond 
20 AU where ionising cosmic rays and X-rays penetrate to 
the mid-plane). Accretion through the “dead zone”, situated 
between these regions of active turbulence, can occur in 


ionised surface layers far above the mid-plane (Oishi et al. 


2007^, from disk winds [Bai and Stone\ 20131 and by purely 


hydrodynamic al instabilities in the vertical shear of the gas 
^Nelson et al.\ |2013| l or radial convection arising from the 
subcritical baroclinic instability {Klahr and Bodenheimer 


2003 Lesur and Papaloizou\ 20101. The mid-plane is be¬ 
lieved to be stirred to a mild degree by these hydrodynam- 
ical instabilities or by perturbations from the active layers 
several scale-heights above the mid-plane, driving effective 
turbulent diffusivities in the interval from 10“^ to 10“^ in 
the mid-plane. The inner and outer edges of this “dead 
zone”, where the turbulent viscosity transitions abruptly, 
are also possible sites of pressure bumps and large-scale 


Rossby vortices which feed off the pressure bumps {Lyra 
[eFHT] [20081 [2009| . 


6.1. Streaming instability 

The low degree of turbulent stirring in the asteroid for¬ 
mation region also facilitates the action of the streaming 
instability, a mechanism where particles take an active role 
in the concentration process {Youdin and Goodman\\2005\ 


\Youdin and Johansen] [2007 [ \Johansen and Youdin\ |2007| l. 

The instability arises from the speed difference between gas 
and solid particles. The gas is slightly pressure-supported in 
the direction pointing away from the star, due to the higher 
temperature and density close to the star, which mimics a 
reduced gravity on the gas. The result is that the gas orbital 
speed is approximately 50 m/s slower than the Keplerian 
speed at any given distance from the star. Solid particles 
are not affected by the global pressure gradient - they would 
move at the Keplerian speed in absence of drag forces, but 
drift radially due to the friction from the slower moving gas. 
The friction exerted from the particles back onto the gas 
leads to an instability whereby a small overdensity of parti¬ 
cles accelerates the gas and diminishes the difference from 
the Keplerian speed. The speed increase in turn reduces 
the local headwind on the dust. This slows down the ra¬ 
dial drift of particles locally, which leads to a run-away pro¬ 
cess where isolated particles drift into the convergence zone 
and the density increases exponentially with time. This pic¬ 


ture is a bit simplified as Youdin and Goodman (20051 and 
Jacquet et al. (|2011|l showed that the streaming instability 


operates only in presence of rotation, i.e. the instability re 
lies on the presence of Coriolis forces. This explains why 
the instability occurs on relatively large scales of the pro¬ 
toplanetary disk where Coriolis forces are important, typi¬ 
cally a fraction of an astronomical unit, and operates most 
efficiently on large particles with frictional coupling times 
around 1/10 of the orbital time-scale (typically dm sizes at 
the location of the asteroid belt). 

Computer simulations which follow the evolution of the 
streaming instability into its non-linear regime show the 
emergence of axisymmetric filaments with typical separa¬ 
tions of 0.2 times the gas scale-height {Yang and Johansm 


20141 and local particle densities reaching several thousand 


times the gas density {Bai and Stone 2010t Johansen et al. 


20121. These high densities trigger the formation of large 


planetesimals (100-1000 km in diameter) by gravitational 
fragmentation of the filaments { Johansen et al.\ 2007| l, al¬ 
though planetesimal sizes decrease to approximately 100 
km for a particle column density comparable to that of the 
solar protoplanetary disk {Johansen et al. \ 2012| l. 

An important question concerning planetesimal forma¬ 
tion through the streaming instability is whether the pro¬ 
cess can operate for particles as small as chondrules in 
the asteroid belt. In Figu re w e show numerical experi¬ 
ments from II Carrera et ai (20151 on the streaming instabil¬ 
ity in particles with sizes down to a fraction of a millimeter. 
The streaming instability requires a threshold particle mass¬ 
loading Z = Ep(Eg, where 27p and Eg are the particle and 
gas column densities, to trigger the formation of overdense 
filaments {Johansen et al. \ |2009b| \Bai and Stone\ |2010|l. 
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Fig. 3.— Space-time plots of particle concentration by streaming instabilities, from Carrera et al. ( 2015| l, with the x-axis 
indicating the radial distance from the centre of the simulation box and the y-axis the time (on the left) and the dust-to-gas 
ratio (on the right). The four columns show particle sizes 0.5 mm (Tf = 0.00117“^), 1.5 mm (rt = 0.00317“^), 5 mm 
(Tf = 0.0117“^) and 1.5 cm (Tf = 0.0317“^). Simulations start with a mean dust-to-gas ratio of Z = 0.01, but gas 
is removed on a time-scale of 30 orbits (1 orbit = 27rl7“^), increasing the dust-to-gas ratio accordingly . While cm-sized 
particles concentrate in overdense filaments already at a modest increase in dust-to-gas ratio Xa Z = 0.015, smaller particles 
require consecutive more gas removal to trigger clumping. 


The simulations in Figure start ?A Z = Zq = 0.01, but 
the particle mass-loading is continuously increased by re¬ 
moving the gas on a time-scale of 30 orbital periods. This 
was done to identify how the critical value of Z depends on 
the particle size. The result is that overdense filaments form 
already at Z = 0.015 for cm-sized particles, while large 
chondrules of mm sizes require Z = 0.04 to trigger fila¬ 
ment formation. Chondrules smaller than mm do not form 
filaments even at Z = 0.08. 

A lowered gas column density may thus be required to 
trigger concentration of chondrule-sized particles by the 
streaming instability. It is possible that the solar proto¬ 
planetary disk had a lower gas density than what is in¬ 
ferred from the current mass of rock and ice in the planets 
(which multiplied by 100 gives the MMSN), if the planet¬ 
forming regions of the nebula were fed by pebbles drifting 
in from larger orbital distances ^Birnstiel et al. |2012| l. In 
this picture the growing planetesimals and planets are fed 
by drifting pebbles, so that the current mass of the plan¬ 
ets was achieved by the integrated capture efficiency of the 
drifting solids; this allows for gas column densities lower 
than in the Minimum Mass Solar Nebula to be consistent 
with the current masses of planets in the Solar System. The 
gas will also be removed by accretion and photoevaporation 
{ Alexander and Armitage\p.QQ6) . The high mass-loading in 
the gas could be obtained through pile-up by radial drift and 


release of refractory grains near the ice line {Sirono 20111. 

Turbulence as weak as a ~ 10“^ is necessary to allow 
the sedimentation of chondrules (with St ~ 10“^) into a 
thin mid-plane layer with scale-height Hp = O.OliTg and 
Pp « pg (see equation]^, the latter being a necessary den¬ 
sity criterion for activating particle pile up by streaming in¬ 
stabilities. Very low levels of a are consistent with proto¬ 
planetary disk models where angular momentum is trans¬ 
ported by disk winds and the mid-plane remains laminar 
2013|l, except for mild stirring by Kelvin- 


{Bai and Stone 


Helmholtz ( Youdin and Shu 2002} and streaming instabili¬ 
ties ^Bai and Stone 2010| l. 

Weak turbulence also facilitates the formation of dm- 
sized chondrule aggregates {Ormel et al. 2008|l, which 


would concentrate much more readily in the gas. Stirring 
by hydrodynamical instabilities in the mid-plane, such as 
the vertical shear instability ^Nelson et al. |2013|), would 


preclude significant sedimentation of chondrule-sized par¬ 
ticles and affect the streaming instability, as well as the for¬ 
mation of chondrule aggregates, negatively. An alternative 
possibility is that the first asteroid seeds in fact did not form 
from chondrules (or chondrule aggregates), but rather from 
larger icy particles which would have been present in the 
asteroid formation region in stages of the protoplanetary 
disk where the ice line was much closer to the star ( |Mar-| 
\tin and Livio\\20\2\^os and Johans^\20\2>}. Chondrules 
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could have been incorporated into by later chondrule accre¬ 
tion (see Section]^. 

7. LAYERED ACCRETION 

The turbulent concentration model and the streaming in¬ 
stability, reviewed in the previous sections, are the leading 
contenders for primary accretion of chondrules into chon¬ 
drites. However, none of the two are completely successful 
in explaining the dominance of chondrules in chondrites: 
the turbulent concentration models may not be able to con¬ 
centrate sufficient amounts for gravitational collapse, while 
the streaming instability relies on the formation of chon¬ 
drule aggregates and/or gas depletion and pile up of solid 
material from the outer parts of the protoplanetary disk. In 
the layered accretion model the chondrules are instead ac¬ 
creted onto the growing asteroids over millions of years af¬ 
ter the formation of the first asteroid seeds - those first seeds 
forming by direct coagulation from a population of 100-m- 
sized planetesimals as envisioned in Weidenschilling ( |201 1 1 
or by one or more of the particle concentration mechanisms 
described in the previous sections. 


7.1. Chondrule accretion 

Chondrules are perfectly sized for drag-force-assisted 
accretion onto young asteroids. The ubiquity of chon¬ 


drules inside chondrites, and their large age spread {Con¬ 


nelly et al. 20121, indicates that planetesimals formed and 
orbited within a sea of chondrules. Chondrules would 
have been swept past these young asteroids with the sub- 
Keplerian gas. The gas is slightly pressure-supported in 
the radial direction and hence moves slower than the Ke- 
plerian speed by the positive amount Aw ^Weidenschillin g 
|1977b[ pVakagawa et al. 19861. The Bondi radius i?B = 
GM/{Av)‘^ marks the impact parameter for gravitational 
scattering of a chondrule by an asteroid of mass M, with 




R 


50 km 


Av 

53m/s 


3.5g/cm3 


(4) 

Here we have normalised by Av = 53 m/s, the nominal 
value in the Minimum Mass Solar Nebula model of\HayashJ\ 


( 1981), and used the chondrule density p, = 3.5 g/cm^ as a 


reference value. Chondrules with friction time comparable 
to the Bondi time-scale /b = Rb/A v are accreted by the 


asteroid ( 

Johansen and Lacerda 

20101 \Ormel and Klahr\ 

\20l0\\Lambrechts and Johansen 

201211. The accretion ra- 


dius Race can be calculated numerically as a function of 
asteroid size and chondrule size by integrating the trajec¬ 
tory of a chondrule moving with the sub-Keplerian gas flow 
past the asteroid. The accretion radius peaks at Race ~ Rb 


for ti/tB in the range from 0.5 to 10 (Lambrechts and Jo 


hansen\\20\2\. Accretion at the full Bondi radius happens 


for particle sizes 


a = [0.008,0.16] mm 


f R \ 


Av 


-3 


\^50kmy \53m/s^ 

yf ^ rV V (5) 

V2.5AU/ V^mmsn/ 

An asteroid of radius 50 km thus “prefers” to accrete chon¬ 
drules of sizes smaller than 0.1 mm, corresponding to the 
smallest chondrules found in chondrites. At 100 km in ra¬ 
dius the preferred chondrule size is closer to 0.2 mm, a 200 
km radius body prefers mm-sized chondrules and larger 
bodies can only grow efficiently if they can accrete chon¬ 
drules of several mm or cm in diameter. Carbonaceous 
chondrites accreted significant amounts of CAIs and ma¬ 
trix together with their chondrules; Rubin ( 2011| l suggested 
that matrix was accreted in the form of cm-sized porous 
aggregates with aerodynamical friction time comparable to 
chondrules and CAIs. 

Aerodynamical accretion of chondrules could explain 
the narrow range of chondrule sizes found in the various 
classes of meteorites. The model predicts that asteroids ac¬ 
crete increasingly larger chondrules as they grow. This pre¬ 
diction may be at odds with the little variation in chondrule 
sizes found within chondrite classes (70% of EH3 and C03 
chondrules have apparent diameters within a factor of 2 of 


the mean apparent diameters in the group, according to Ru- 
2000|l. The least metamorphosed LL chondrites nev- 


bin 


ertheless do seem to host on the average larger chondrules 
{Nelson and Rubin 2002 | l. More metamorphosed LL chon¬ 
drites actually show a lack of small chondrules; this could 
be due to the fact that the smallest chondrules disappeared 
from the strongly heated central regions of the parent body. 


The accretion rate of chondrules (and other macroscopic 
particles) is 

M = TTf^R^ppAv. ( 6 ) 

Here /b parameterises the actual accretion radius relative to 
the Bondi radius and pp is the chondrule density. Accretion 
of chondrules is a run-away process, since M cx i?| cx 
if the optimal chondrule size is present (so that /b = 1 in 
equation]^. The characteristic growth time-scale is 


M 

/exp = ^ = 1.66Myr 

M 


R 


50 km 
-1 


-3 


Av 

53m/s 
-1 


V (7) 

V2.5AUy' \SmmsnJ \^.5g/cm^J ^ ^ 


We assumed here that the chondrules have sedimented to 
a thin mid-plane layer of thickness 1% relative to the gas 
scale-height. The strong dependence of the accretion rate 
on the planetesimal mass will drive a steep differential 
size distribution of a population of planetesimals accreting 
chondrules. This is illustrated in Figure where asteroid 
seeds with initial sizes from 10 to 50 km in radius have been 
exposed to chondrule accretion over 5 Myr. The value of the 
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Fig. 4.— The size distribution of asteroids and planetary 
embryos after accreting chondrules with sizes from 0.1 to 
1.6 mm in diameter for 5 Myr. The original asteroid sizes 
had sizes between 10 and 50 km in radius (yellow line), 
here envisioned to form by the streaming instability in a 
population of dm-sized icy particles. The resulting size dis¬ 
tribution of asteroids is in good agreement with the bump at 
70 km in radius, the steep size distribution from 70 km to 
200 km, and the shallower size distribution of larger aster¬ 
oids whose chondrule accretion is slowed down by friction 
within their very large Bondi radius. Figure based on the 
results Qf \Johansen et q/.| ( |2015| |. 

turbulent viscosity is a = 2 x 10“®. The resulting differen¬ 
tial size distribution (which manifests itself after around 3 
Myr of chondrule accretion) shows a bump at 70 km in ra¬ 
dius, a steep decline towards 200 km in radius, and finally a 
more slow decline towards larger asteroids. The shallower 
decline is caused by the lack of cm-sized particles needed 
to drive the continued run-away accretion of large asteroids 
(this can be seen as a drop in /b in equation]^. All the fea¬ 
tures in the size distribution in Figure are in good agree¬ 
ment with features of the observed size distribution of as¬ 
teroids which are not explained well in coagulation models 
^Morbidelli et a/.]|2009| ). 

Layered accretion of chondrules can readily explain the 
large age spread of individual chondrules inside chondrites 
( Connelly et a/.[ |2012|l, as well as the remanent magnetisa¬ 
tion of the Allende meteorite \Elkins-Tanton et ar\ |201 l| l, 
imposed on the accreted chondrules from the internal dy¬ 
namo in the parent body’s molten core. Efficient chon¬ 
drule accretion requires, as does the streaming instability 
discussed in the previous section, sedimentation of chon¬ 
drules to a thin mid-plane layer. The turbulent viscosity of 
a = 2 X 10“® in Figure|^is nevertheless significantly larger 
than the a ~ 10”^ needed to sediment chondrules to a thin 


mid-plane layer of thickness 1% of the gas scale height; and 
indeed the time-scale to grow to the current asteroid popu¬ 
lation is 3 Myr in the simulation shown in Figure 4, about 
twice as long as in equation Q, but in good agreement with 
the ages of the youngest chondrules. 

8. OPEN QUESTIONS 

The formation of asteroids is a complex problem which 
will only be solved through a collective effort from as¬ 
tronomers, planetary scientists and cosmochemists. Al¬ 
though many details of asteroid formation are still not un¬ 
derstood, we hope to have convinced the reader that new 
insights have been achieved in many areas in the past years. 
Flere we highlight ten areas of open questions in which we 
believe that major progress will be made in the next decade: 

1) Short-lived radionuclides. What is the origin of the 
short-lived radioactive elements which melted the differen¬ 
tiated parent bodies? Was ^®A1 heterogeneous in the solar 
protoplanetary disk? How did the young Solar System be¬ 
come polluted in ^®A1 without receiving large amounts of 
eope, 

an element that is copiously produced in supernovae? 

2) Maintaining free-floating chondrules and CAIs. How 
is it possible to preserve chondrules and CAIs for millions 
of years in the disk before storing them in a chondritic body, 
without mixing them too much to erase chondrule classes 
and chondrule-matrix complementarity? What are we miss¬ 
ing that makes this issue so paradoxical? 

3) Chondrules versus matrix. Why do carbonaceous 
chondrites contain large amounts of matrix while ordinary 
chondrites contain very little matrix? Did the matrix en¬ 
ter the chondrites as (potentially icy) “matrix lumps” or on 
fine-grained rims attached to chondrules and other macro¬ 
scopic particles? 

4) Initial asteroid sizes. What is the origin of the steep 
differential size distribution of asteroids beyond the knee 
at 100 km? Did asteroids form small as in the coagulation 
picture, medium-sized as in the layered accretion model or 
large as in some turbulent concentration models? Why do 
Kuiper belt objects, which formed under very different con¬ 
ditions in temperature and density, display a similar size 
distribution as asteroids? 

5) The origin of asteroid classes. How is the radial gra¬ 
dient of asteroid composition produced and retained in the 
presence of considerable pre-accretionary turbulent mixing 
and post-accretionary dynamical mixing? Is asteroid for¬ 
mation a continuous process which happens throughout the 
life-time of the protoplanetary disk? What do the different 
chondrite groups mean in terms of formation location and 
time? 

6) Dry and wet chondrites. Why do we have dry chon¬ 
drites (enstatite, ordinary)? If chondrites form at 2-4 Myr 
after CAIs, then the snowline should have been well in¬ 
side the inner edge of the asteroid belt. Are there over¬ 
looked heating sources which could keep the ice line at 3 
AU throughout the life-time of the protoplanetary disk? Or 
did the asteroid classes form at totally different places only 
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to be transported to their current orbits later? 

7) Internal structure of asteroids. Does the chondrite and 
asteroid family evidence suggest that the primary asteroids 
- before internal heating - are homogeneous, roughly 100- 
km-diameter bodies composed of a physically, chemically, 
and isotopically homogeneous mix of chondrule-sized com¬ 
ponents? Or is internal heterogeneity, as may be the case for 
the Allende parent body, prevalent? 

8) Turbulent concentration of chondrules. Under what 
nebula conditions can vortex tubes over a range of nebula 
scales concentrate enough chondrules into volumes that are 
gravitationally bound, at a high enough rate to produce the 
primordial asteroids and meteorite parent bodies directly? 
What other roles could turbulent concentration play in plan- 
etesimal formation given that the optimally concentrated 
particle is chondrule-sized under nominal values of the tur¬ 
bulent viscosity? 

9) Streaming instability with chondrules. Will the con¬ 
ditions for streaming instabilities to concentrate chondrule- 
sized particles, i.e. gas depletion and/or particle pile up, be 
fulfilled in the protoplanetary disk? How does an over- 
dense filament of chondrule-sized particles collapse under 
self-gravity given the strong support by gas pressure? 

10) Layered accretion. What is the origin of the apparent 
scarcity of heterogeneous asteroid families, given that aster¬ 
oids orbiting within an ocean of chondrules should accrete 
these prodigiously? What is the thermal evolution of early- 
formed asteroid seeds that continue to accrete chondrules 
over millions of years? 
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